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Abstract 



We investigate the thermodynamic arrow of time in a time-symmetrically 
recollapsing universe by calculating quantum mechanically the entropy pro- 
duction of a massive scalar field. It is found that even though the Hamiltonian 
has a time-reversal symmetry with respect to the maximum expansion of the 
universe, the entropy production is generic and the total entropy of the scalar 
field increases monotonically. We conclude that the thermodynamic arrow 
of time is a universal phenomenon even in the expanding and subsequently 
recollapsing universe due to the parametric interaction of matter field with 
gravity. 
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I. INTRODUCTION 



The thermodynamic arrow of time has always been a puzzling issue in an expanding and 
subsequently recollapsing universe. When the Universe is to recollapse time-symmetrically 
with respect to its maximum expansion, can the entropy of matters always increase or time- 
symmetrically decrease? Historically, Gold gave the first scientific explanation that the 
entropy of matter should decrease time-symmetrically with respect to the maximum expan- 
sion of the universe, but this argument was refuted by Penrose 0. In the context of quantum 
cosmology, Hawking || argued that the entropy of matter should be time symmetric since 
the Hartle and Hawking's no-boundary wave function || should be CTP invariant. But Page 
H refuted this argument by pointing out the time-asymmetric wave function obeying CTP 
theorem, in which the entropy can increase monotonically throughout the expansion and 
subsequent recollapse of the universe. Later Hawking et al. || admitted the thermodynamic 
arrow of time by considering the growing density perturbation even through the recollapsing 
universe. Recently, we have argued that the parametric interaction between the gravity and 
matter fields may break not only the cosmological symmetry of time @ but also the ther- 
modynamic symmetry of time M . Kiefer and Zeh || also argued the thermodynamic arrow 
of time in the recollapsing universe. (See Ref. ]10| for comprehensive review and references.) 

In this paper we revisit the issue of thermodynamic arrow of time by studying a massive 
scalar field in a time-symmetrically recollapsing universe. Matter field in such a time- 
dependent background spacetime should be described by a nonequilibrium quantum field. 
We apply the recently introduced unified approach [[Lip based on the Liouville-von Neumann 
(LvN) equation to the nonequilibrium massive scalar field and find the time-dependent Fock 
(Hilbert) spaces throughout the expansion and subsequent recollapse of the universe. The 
particle production is calculated explicitly and the information-theoretic entropy is computed 
for the nonequilibrium quantum process. From this we conclude that the thermodynamic 
arrow of time is a universal phenomenon due to the parametric interaction of the matter 
field with the gravity even in the recollapsing universe. 

The organization of this paper is as follows. In Sec. II we review the nonequilibrium 
quantum evolution of the massive scalar field in the recollapsing Friedmann-Robertson- 
Walker (FRW) universe. In Sec. Ill we study an exactly solvable model for the massive 
scalar field in the time-symmetrically recollapsing universe. Finally, in Sec. IV we study the 
thermodynamic arrow of time in the recollapsing FRW universe. 



II. MASSIVE SCALAR FIELD IN THE FRW UNIVERSE 



In order to study the thermodynamic arrow of time in the time-symmetrically recollaps- 
ing universe, we consider the simple model of a closed FRW universe minimally coupled to 
a massive scalar field. The massive scalar field is described by the action 



/ = — - / (l'.V\'- 



The FRW universe has the metric 

ds 2 = -dt 2 + R 2 {t)dQ 2 3 , 



(1) 



(2) 
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where d£l\ is the metric on the unit three-sphere. According to Ref. [12] the scalar field can 
be decomposed as 



S(t,x)= E#Jf)Q? m (x), ( 3 ) 



n,Lm 



where Q" m denote the eigenf unctions of the Laplace operator on the three-sphere and are 
explicitly given by Qf m = n™Y; m , where Y^ m are the spherical harmonics on the two-sphere 
and IT™ are the Fock harmonics. Then from the action ([[J) we obtain the mode-decomposed 
Hamiltonian for the scalar field 

where k denote collectively in, I, m) and tt^ are momenta conjugate to (j) k . Thus, the Hamil- 
tonian is a collection of infinitely many time-dependent harmonic oscillators. 

Each time-dependent harmonic oscillator in Eq. ([|) describes a nonequilibrium quan- 
tum process through the parametric interaction between the matter and gravity. Such 
time-dependent system leads necessarily to the particle production in quantum mechanics 
and its nonequilibrium process implies the entropy production in quantum statistical me- 
chanics. To treat the nonequilibrium process appropriately, we shall adopt the recently 
introduced unified approach (JXTJ] based on the LvN equation. This LvN approach provides 



us explicitly not only with the exact quantum states but also with the density operator. In 
this LvN approach the fundamental laws for the nonequilibrium quantum system are the 
time-dependent Schrodinger equation (in unit h = 1) 

and the LvN equation for the density operator 



dt 



+ [p(t),H(t)]=0. (6) 



It was Lewis and Riesenfeld ||T3[ who observed that the eigenstate of any operator O satisfy- 
ing the LvN equation (^) is an exact quantum state of the time-dependent quantum system 
up to a trivial time-dependent phase factor. (For development and various applications, see 
Ref. @.) 



The stratagem of the LvN approach is to find the pair of time- dependent operators 

A k (t) = i[u* k (t)n k - R 3 (t)u k (t)4> k ], AUt) = h.c.(i k (t)), (7) 
that satisfy the LvN equation 

.dAv(t) 



i- 



dt 



+ [A lt (t),H lt (t)] = 0. (8) 



It readily follows that Eq. (H) is satisfied only when u k is a complex solution to the corre- 
sponding classical equation of motion 
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R(t), 



rr 



u k {t) = 0. 



(9) 



A k (t) and A£,(t) can be made the annihilation and creation operators with the standard 
commutation relation at each time t 



i. 



[A k {t),Al{t)\ =8 Kkl 

by imposing the Wronskian condition 

R 3 (t) [u k {t)ut{t) - K(t)u k (t) 
Then, neglecting trivial time- dependent phase factors, the number states 

n k ,t = -=(A k (t ) k ,t), 

together with the vacuum state 

A k (t)\0 k ,t) = 0, 



(10) 



(11) 



(12) 



(13) 



constitute the time-dependent Fock (Hilbert) space. In fact, there are one-parameter Gaus- 
sian states flHf ) for each oscillator and the vacuum state is the Gaussian state having the 
minimum uncertainty [13]. 

Note that the Hamiltonian system (|j) may not be solved exactly since the classical 
equations of motion (|S[) can not be solved in closed form for the scale factor 



R(t) 



3m 2 p 

csch 



\ 8np 



(14) 



where nip = ^= is the Planck mass and p is the density of dust particles in the the closed 
FRW universe. So it would be interesting to study first a model that can be solved exactly. 



III. AN EXACTLY SOLVABLE MODEL 

We now consider a model Hamiltonian that shares many physically interesting features 
with Eq. (f|). The frequencies of oscillators in Eq. @ increase indefinitely as the universe 
(re)collapses (R — > 0), and approach to constant values at the maximum expansion. To 
study the particle production in a curved spacetime, the following Hamiltonian has been 
often used [[L5[ 



m 0f,,2 , ^ 2 N 



2mo 

whose the classical equations of motion can be solved analytically 

,2 



(15) 



^) + K + ^>k(*) 



0. 



(16) 
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In terms of the complex solution for each oscillator one finds the annihilation and creation 
operators 



B k {t) = * [<(t)Pk ~ mot£(*)fflc] , B{(t) = h.c. (S k (t)) . 
Further, one should impose the Wronskian condition 



to 



v k (t)vt(t) - v* k (t)u k (t) 



(17) 



(18) 



to guarantee the standard commutation relation [I? k (t), B^,(t)] = 5 kjk / for all times. 
The complex solution to Eq. (|TE|) satisfying the Wronskian (113) is given by 



(19) 



where W(a, ±r) are the parabolic cylinder functions and 



a k = 9-, r = k = \fl + e 2 ™" _ e ™k_ 

In the past infinity (t — > — oo) the solution ( [Of) has the asymptotic form [16| 



exp 



T -a k ln(-r) + I + y 



(20) 



(21) 



where tp% = argr(J + ia k J. Note that the asymptotic solution fl2"l|) will be the same as the 
adiabatic (WKB) solution of Sec. IV. Also another useful form of the solution ( |19"D is given 
by 



k--)E(a k ,T) + l -(k + ^)E\a^r 



(4mgu*)V4 

In the future infinity (t — > oo), u kout has the asymptotic form 

1 



(22) 



(4m^u; k ) 1 /4 



(4mgu; k T 



2^1 1/4 



exp 



—i 



[— - a k ln(r) + - + 



7T (£2 

4 2 



(23) 



From Eqs. (|2l| ) and (|23"D we are able to find the transformation of solutions between two 
asymptotic regimes 



Vk, in (t) = ~(k + -)v k>out {t) + -(«- ^)v k)OUt (t), 
and the Bogoliubov transformation 

B ki out{t) = HkBk, in (t) + is k Bl in (t), h.c. 

where 



(24) 



(25) 
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Also the Bogoliubov transformation can be rewritten as 

B^ out {t) = S k (z k )B Kin (t)Sl(z k ), h.c. 



7T Af 

k,in- D k,in 



exp 



rk m 2 

2 v-^k^n 



in terms of the squeeze operator [I7| 

Sk(>k) = exp 
in which 

/Uk = coshrke'a, z/ k = sinhr k e*~ 



Rt2 s 
■°k,in; 



(26) 



(27) 



(28) 



(29) 



The vacuum states (fL3|) for each oscillator depend on time explicitly and those at two 
different times are related by the Bogoliubov transformation 



where 



B k (t) = /i k (t,t )B k (t ) + u k (t,t )Bl(t ), h.c 



jUk(Mo) = im [vl(t)v k {to) - vl(t)v k (t )], 
Vk{t,t ) = im [vl(t)vl(t ) - vi(t)vl(t )}. 



(30) 



(31) 



The quantum state of each oscillator evolves unitarily by the unitary transformation (p7|). 
and the vacuum states at two different times are not orthogonal to each other 



(Ok,t|O k ,*o) 



A(t>*o) \ 1/2 



(wmo)|, 



(32) 



where i?k = arg(f k (^o)' y k(^))- However, since fi^ > 1 and can be unity at most for a finite 
number of modes, the vacuum states for all modes of the scalar field 



io,t> = n io k ,t>, 



(33) 



are orthogonal to each other at two different times 

e itf k (Mo) \ 1 / 2 



(0,t|0,t ) = n 

k 



(|/*k(Mo) 



(34) 



Therefore, the time-dependent Fock spaces constructed from the vacuum states form uni- 
tary inequivalent representations. The physical implication of the unitary inequivalence of 
representations is that particles are produced in the infinite future (t — > oo) from the initial 
vacuum state of the infinite past (t — > — oo) by the amount 



H B i,out B *,out\®)in = u l- 
k k 



(35) 
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Similarly, the particle production from the number states is found to be 
[IL«K|] Bi,outBk,out[Y[\n k ) in ] =XX( 2 "k+ 1). 



(36) 



The particle production in the expanding universe was first explained by Parker [T|| . 

We can make use of the operators flTTj ) satisfying the LvN equation to construct a density 
operator 



Pk(t) 



1 



Z k (t) 



-^n k (^(t)J3 k (t)+|) 



Z k (t) = Tre 



-^n k (B k (t)fl k (t)+i) 



(37) 



Here /3 is a free parameter for nonequilibrium quantum processes, which becomes (3 = 
for an equilibrium system, where k B and T are the Boltzmann constant and temperature, 
respectively. Using the information-theoretic entropy at each time defined by |T9| 



Sk,ent.{t) = -k B Ti p k (t) ln/S k (t) 



(38) 



the entropy production during the evolution from the infinite past to the infinite future is 
given by 



k,ent. 



k B {Tr 



Pk,out hi p\ 



out 



Tr 



P\,in hi p\ ■ 



(39) 



It should be pointed out that the entropy cannot be produced for a system of finite oscillators 



due to the unitary transformation (27). 

However, the Fock space of the scalar field in the infinite future is unitarily inequivalent 
to that of the infinite past. In other words, Eq. (|27D should not be regarded as a unitary 
transformation but rather as a unitarily inequivalent transformation of the scalar field. 
Therefore, the entropy production for the scalar field should be evaluated with respect 
to the same Fock space, for instance, that of the infinite past 



AS] 



U.,ent. 



Pout In p out \n k ) 

in in Pin 

In p in |n k ) 

in ( 

-k B I m(^k|Anlnpi„|ri k )i n (|i n (rn k |5'|n k )i n | 2 - 8. 



m k ,n k 



(40) 



n k ,m k 



Here we made use of the transformation (|27j ) in the second line. For a small squeeze param- 
eter r k we obtain approximately AS^ ent . — 2r k 5 , kj j„, which is different from AiS , kjenfc ~ 2r k 
for the Gaussian states in Refs. pO| , pT| . The total amount of the information-theoretical 
entropy is given by 



AS ent . ~ £ 2 hi(V 1 + e~ 2 ^ - e"^)^. 



(41) 



It is remarkable that the entropy increases monotonically during the symmetric evolution 
from the infinite past to the future and the amount of the entropy production is proportional 
to the entropy present in the past. 
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IV. THERMODYNAMIC ARROW OF TIME IN THE RECOLLAPSING 

UNIVERSE 



We now return to the scalar field in the FRW universe and its Hamiltonian (|]). Oft- 
times the instantaneous annihilation and creation operators are used that diagonalize the 
Hamiltonian at every moment: 



A ktdig .(t) = e'/<W*> [ ^ 3( ^ k(t) g k + ^ 3( ^ k(t) P k ] , iU,W = h.c.(i Mis .(*))- (42) 

However, the Fock spaces of these number states lead to the infinite production of particle 
and entropy On the other hand, the LvN approach in Sees. II and III will yield the 
correct and physically meaningful result, provided that the analytical solution of Eq. @ 
is used for each mode. Since it is difficult to solve analytically Eq. (|9J) for the scale factor 
(|14T), we rely on the adiabatic (WKB) solution. This is done by interpreting Eq. (^J) as a 
one-dimensional Schrodinger equation and, for a slowly varying frequency u>k(t), by using 
the adiabatic solution of the form 



u k ,ad(t) = =exp -i / Q k (t) . (43) 



The adiabatic annihilation and creation operators for each oscillator are related with the 
instantaneous ones via the Bogoliubov transformation 



^-W = i 1 - - isr^'^UM, - (44) 



As for the model in Sec. Ill, the adiabatic solution of the infinite past can be written in 
terms of that of the infinite future as 

The explicit forms of /ik and depend on R(t), m, and k. Similarly, the adiabatic solution 
at t p in the expanding stage can be expressed in terms of that at tf in the time-symmetric 
recollapsing stage as 

u^ad.{t P ) = ^k(t f ,t p )u Kad Xtf) - Vk(tf,t p )u1 ad (tf). (46) 

In general, one has fik(tf,t p ) ^ even for symmetric times t p and tf with respect to the 
maximum expansion of the universe. One can also understand this fact from the scattering 
of a particle by a symmetric potential in quantum mechanics. From the study of the model 
in Sec. Ill we may find some characteristic features of Eq. (f|). It is rather straightforward 
to obtain the entropy production for each mode of the scalar field 

&S Kent (t f ,t p ) ~ 2r k (t f ,t p )S k (t p ), (47) 

where 

cosh r k (t f , t p ) = \fi k (t f ,t p )\. (48) 
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The total entropy production of the scalar field is then given by 

AS ent Xt f ,t p ) -J2 2r ^( t f^p) S M- ( 49 ) 

k 

Thus we have shown that the entropy of the matter (scalar) field increases even during the 
time-symmetric recollapse of the universe. This thermodynamic arrow of time is due to 
the parametric interaction between the matter field and gravity. The particle and entropy 
production seems to be a generic feature and thereby, the global thermodynamic arrow of 
time is a consequence of the evolution of the universe. 



V. CONCLUSION 



In this paper we have studied the particle and entropy production of a massive scalar 
field in a time-symmetrically recollapsing universe. In the expanding and subsequently rec- 
ollapsing universe the matter field is described by a nonequilibrium quantum field due to the 
parametric interaction of the matter field with the gravity. To treat properly we have applied 
the unified approach based on the Liouville-von Neumann equation to the nonequilibrium 
massive scalar field and found the time-dependent Fock spaces throughout the evolution of 
the universe. It is the parametric interaction that leads to the entropy production during the 
evolution of the universe. The study of an exactly solvable model which has many similar 
features with the scalar field revealed the monotonic increase of particle and entropy even 
for the time-symmetric recollapse of the universe. Likewise, the entropy of a massive scalar 
field in the expanding and recollapsing universe is found to be monotonically increasing 
using adiabatic (WKB) solution to the classical equation of motion. Therefore, we conclude 
that the entropy of the scalar field increases monotonically throughout the evolution of the 
universe including even the symmetrically recollapsing universe. 

In this paper we have focused only on the scalar field in the expanding and subsequently 
recollapsing universe. The other matter of the Universe is the fermionic field. Recently, the 
Liouville-von Neumann approach has been applied even to the time-dependent fermionic 
field to describe correctly its nonequilibrium evolution P^| . Since the coupling constants of 
the fermionic field depend on time as the Universe expands and subsequently recollapses, 
the result of Ref. J22| can be employed to study the thermodynamic arrow of time due to 
the fermionic field, which will be presented in a future publication. 
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